**************************************************************************************************

******** The effect of COVID-19 lockdowns on political support: Some good news for democracy? ****


* Damien Bol, Marco Giani, André Blais, Peter John Loewen


                *REPLICATION OF MAIN FIGURE* 
				
	
ssc install ebalance
ssc install coefplot
ssc install estout

clear all
use Data.dta


**STFdem
 
quietly: reg Stfdem2 Treatment  i.Country, vce(cluster  CL_var) 
est store Basic
quietly: reg Stfdem2 Treatment  Deaths  Date    i.Country, vce(cluster  CL_var)
est store Augmented
quietly: reg Stfdem2 Treatment  Deaths  Date Age Female Rural Immigrant University Vote i.Country, vce(cluster CL_var)
est store Full
ebalance Treatment  Age Female Rural Immigrant University Vote 
quietly: reg Stfdem2 Treatment  Deaths  Date Age Female Rural Immigrant University Vote  i.Country [pweight=_webal], vce(cluster CL_var)
est store Fully


	   grstyle init
	 grstyle set legend 2,  nobox
	 grstyle set size 8pt: tick_label key_label
	 grstyle set size 12pt: heading
	 grstyle set size 10pt: subheading axis_title
	grstyle set color   "224 224 224" "178 178 178" "100 100 100"  "0 0 0" 
    grstyle set graphsize 14.5cm 13cm
	coefplot  (Basic, label("Basic" "model")) (Augmented, label("Expanded" "model"))  (Full, label("Full" "model")) (Fully, label("Balanced" "model")) ,  drop(  _cons   *.Country )       xline(0,  lpattern(dash)   ) msymbol(s)  levels(95) ciopts(recast(. rcap))      ylabel(, grid)   headings( Treatment="{bf:Policy change}"  Deaths="{bf:COVID-19 incidence}" Age ="{bf:Socio-demographics}" Vote ="{bf:Political behavior}", labgap(-90) )   legend(position(12) rows(1))    yscale(alt noline) graphregion(margin(l=34))   coeflabels(, labgap(-85) notick)    baselevels   note("Effective sample: N=6,348 (Full model).", span ) 
 

**Trust

quietly: reg Trust2 Treatment  i.Country, vce(cluster CL_var)
est store Basic
quietly: reg Trust2 Treatment  Deaths  Date i.Country, vce(cluster CL_var)
est store Augmented
quietly: reg Trust2  Treatment  Deaths  Date  Age Female Rural Immigrant University Vote  i.Country, vce(cluster CL_var)
est store Full
quietly: reg Trust2 Treatment  Deaths  Date Age Female Rural Immigrant University Vote  i.Country [pweight=_webal], vce(cluster CL_var)
est store Fully


	   grstyle init
	 grstyle set legend 2,  nobox
	 grstyle set size 8pt: tick_label key_label
	 grstyle set size 12pt: heading
	 grstyle set size 10pt: subheading axis_title
	grstyle set color  "224 224 224" "178 178 178" "100 100 100"  "0 0 0"  
	 graph set window fontface default
	 grstyle set graphsize 14.5cm 13cm
	coefplot  (Basic, label("Basic" "model")) (Augmented, label("Expanded" "model"))  (Full, label("Full" "model")) (Fully, label("Balanced" "model")) ,  drop(  _cons   *.Country )       xline(0,  lpattern(dash)   ) msymbol(s)  levels(95) ciopts(recast(. rcap))      ylabel(, grid)   headings( Treatment="{bf:Policy change}"  Deaths="{bf:COVID-19 incidence}" Age ="{bf:Socio-demographics}" Vote ="{bf:Political behavior}", labgap(-90) )   legend(position(12) rows(1))    yscale(alt noline) graphregion(margin(l=34))   coeflabels(, labgap(-85) notick)    baselevels  note("Effective sample: N=6,275 (Full model).", span ) 
 

**Party

quietly: reg partygov Treatment  i.Country, vce(cluster CL_var)
est store Basic
quietly: reg partygov Treatment  Deaths  Date i.Country, vce(cluster CL_var)
est store Augmented
quietly: reg partygov  Treatment  Deaths  Date  Age Female Rural Immigrant University Vote  i.Country, vce(cluster CL_var)
est store Full
quietly: reg partygov Treatment  Deaths  Date Age Female Rural Immigrant University Vote  i.Country [pweight=_webal], vce(cluster CL_var)
est store Fully


	   grstyle init
	 grstyle set legend 2,  nobox
	 grstyle set size 8pt: tick_label key_label
	 grstyle set size 12pt: heading
	 grstyle set size 10pt: subheading axis_title
	grstyle set color  "224 224 224" "178 178 178" "100 100 100"  "0 0 0"  
	 graph set window fontface default
	 grstyle set graphsize 14.5cm 13cm
	coefplot  (Basic, label("Basic" "model")) (Augmented, label("Expanded" "model"))  (Full, label("Full" "model")) (Fully, label("Balanced" "model")) ,  drop(  _cons   *.Country )       xline(0,  lpattern(dash)   ) msymbol(s)  levels(95) ciopts(recast(. rcap))      ylabel(, grid)   headings( Treatment="{bf:Policy change}"  Deaths="{bf:COVID-19 incidence}" Age ="{bf:Socio-demographics}" Vote ="{bf:Political behavior}", labgap(-90) )   legend(position(12) rows(1))    yscale(alt noline) graphregion(margin(l=34))   coeflabels(, labgap(-85) notick)    baselevels note("Effective sample: N=4,367 (Full model).", span ) 



